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According to the dynamical quasiparticle model (DQPM) - matched to reproduce lattice QCD 
results in thermodynamic limit, - the constituents of the strongly interacting quark-gluon plasma 
(sQGP) are massive and off-shell quasi-particles (quarks and gluons) with broad spectral functions. 
In order to address the electromagnetic radiation of the sQGP, we derive off-shell cross sections of 
qq — > 7* , qq — > 7* + g and qg — > 7*g(g<? — ¥ J*q) reactions taking into account the effective propaga- 
tors for quarks and gluons from the DQPM. Dilepton production in In+In collisions at 158 AGeV is 
studied by implementing these processes into the parton-hadron-string dynamics (PHSD) transport 
approach. The microscopic PHSD transport approach describes the full evolution of the heavy-ion 
collision: from the dynamics of quasi-particles in the sQGP phase (when the local energy density is 
above ~ 1 GeV/fm 3 ) through hadronization and to the following hadron interactions and off-shell 
propagation after the hadronization. A comparison to the data of the NA60 Collaboration shows 
that the low mass dilepton spectra are well described by including a collisional broadening of vector 
mesons, while the spectra in the intermediate mass range are dominated by off-shell quark-antiquark 
annihilation, quark Bremsstrahlung and gluon-Compton scattering in the nonperturbative QGP. In 
particular, the observed softening of the mj spectra at intermediate masses (1 GeV < M < 3 GeV) 
is approximately reproduced. 

PACS numbers: 25.75.Nq, 24.10.Lx, 25.75.Cj, 14.40. Be, 24.85. +p 

I. INTRODUCTION 

While the properties of hadrons are rather well known in free space (embedded in a nonperturbative QCD vacuum) 
the masses and lifetimes of hadrons in a baryonic and/or mesonic environment are subject of current research in 
order to achieve a better understanding of the strong interaction and the nature of confinement. A modification 
of vector mesons has been seen experimentally in the enhanced production of lepton pairs above known sources in 
nucleus- nucleus collisions at Super-Proton-Synchroton (SPS) energies [HQ- This can be attributed to a shortening 
of the lifetime of the vector mesons p, uj and <fr. The question arises if the enhancement might (in part) be due to 
new radiative channels Q from the strong Quark-Gluon Plasma (sQGP). The answer to this question is nontrivial 
due to the nonequilibrium nature of the heavy-ion reactions and covariant transport models have to be incorporated 
to disentangle the various sources that contribute to the final dilepton spectra seen experimentally. 

II. THE PHSD APPROACH 

To address the vector meson properties in a hot and dense medium - as created in heavy-ion collisions - we employ 
an up-to-date relativistic transport model, i.e. the Parton Hadron String Dynamics [J] (PHSD) that incorporates 
the relevant off-shell dynamics of the vector mesons as well as the explicit partonic phase in the early hot and dense 
reaction region. PHSD consistently describes the full evolution of a relativistic heavy-ion collision from the initial hard 
scatterings and string formation through the dynamical dcconfincmcnt phase transition to the quark-gluon plasma 
(QGP) as well as hadronization and to the subsequent interactions in the hadronic phase. 

In the hadronic sector PHSD is equivalent to the Hadron-String-Dynamics (HSD) transport approach that has 
been used for the description of pA and AA collisions from SIS to RHIC energies and has lead to a fair reproduction 
of hadron abundances, rapidity distributions and transverse momentum spectra. In particular, HSD incorporates 
off-shell dynamics for vector mesons - according to Refs. § - and a set of vector-meson spectral functions 
that covers possible scenarios for their in-medium modification. The transport theoretical description of quarks and 
gluons in PHSD is based on a dynamical quasiparticle model for partons matched to reproduce lattice QCD results 
in thermodynamic equilibrium (DQPM). The transition from partonic to hadronic degrees of freedom is described by 
covariant transition rates for the fusion of quark-antiquark pairs to mesonic resonanses or three quarks (antiquarks) 
to baryonic states. 

Various models predict that hadrons change in the (hot and dense) nuclear medium; in particular, a broadening 
of the spectral function or a mass shift of the vector mesons has been expected. Furthermore, QCD sum rules 
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FIG. 1: Diagrams contributing to the dilepton production from a QGP: (a) Drell-Yan mechanism, (b) gluon-Compton scattering 
(GCS), (c) vertex correction, (d) gluon Bremsstrahlung (NLODY), where virtual photons (wavy lines) split into lepton pairs, 
spiral lines denote giuons, arrows denote quarks. In each diagram the time runs from left to right. 

indicated that a mass shift may lead to a broadening and vice versa Q; therefore both modifications should be 
studied simultaneously. Thus we explore three possible scenarios: (1) a broadening of the p spectral function, (2) a 
mass shift, and (3) a broadening plus a mass shift. The HSD (PHSD) off-shell transport approach allows to investigate 
in a consistent way the different scenarios for the modification of vector mesons in a hot and dense medium. In the 
off-shell transport, hadron spectral functions change dynamically during the propagation through the medium and 
evolve towards the on-shell spectral function in the vacuum. 

As demonstrated in Ref. [10( the off-shell dynamics is particularly important for resonances with a rather long 
lifetime in vacuum but strongly decreasing lifetime in the nuclear medium (especially u> and </> mesons) but also proves 
vital for the correct description of dilepton decays of p mesons with masses close to the two pion decay threshold. 
For a detailed description of the various hadronic channels included for dilepton production as well as the off-shell 
dynamics we refer the reader to Refs. [ToL Hl|. 



III. DILEPTON RADIATION FROM THE QUARK GLUON PLASMA 



By em ploy ing the HSD approach to the low mass dilepton production in relativistic heavy-ion collisions, it was 
shown in [11] that the NA60 Collaboration data for the invariant mass spectra for p + p~ pairs from In+In collisions 
at 158 A-GeV favored the 'melting p' scenario [14|. Also the data from the CERES Collaboration [15j showed a 
preference for the 'melting p' picture. On the other hand, the dilepton spectrum from In+In collisions at 158 A-GeV 
for M > 1 GeV could not be accounted for by the known hadronic sources (see Fig. 2 o f (ill ). Also, hadronic models 
do not reproduce the softening of the mx distribution of dileptons at M > 1 GeV [lj]. This observation points 
towards a partonic origin but needs further examination. 

The transport theoretical description of the partonic dynamics in PHSD is based on the Dynamical QuasiParticle 
Model (DQPM) [l3| which describes QCD properties in terms of single-particle Green's functions (in the sense of 
a two-particle irreducible approach) and leads to effective strongly interacting partonic quasiparticles with broad 
spectral functions as degrees of freedom. Dilepton radiation by the constituents of the strongly interacting QGP, i.e. 
by the dynamical quasiparticles, proceeds via the elementary processes illustrated in Fig. [T] the basic Drell-Yan q + q 
annihilation mechanism, Gluon Compton scattering (q + g —> 7* + q and q + g — > 7* + q), and quark + anti-quark 
annihilation with gluon Bremsstrahlung in the final state {q + q — > 5 + 7*)- Note that in our calculations the running 
coupling as depends on the local energy density e according to the DPQM [l3j |. while the energy density is related 
to the temperature T by the 1QCD equation of state. Numerically, we observe that as is of the order O(l) and thus 
the contribution of the higher order diagrams in Fig. Q] is not subleading! 

In the on-shell approximation, one uses perturbative QCD cross sections for the processes listed above. However, 
to make a quantitative comparison to experimental data at the realistically low photon virtuality M ~ 1 GeV, one 
has to take into account the non-perturbative spectral functions and self-energies of the quarks, anti-quark and giuons 
thus going beyond the leading twist. Therefore, we have derived the off-shell cross sections of qq — > 7*, qq — > 7* + g 
and qg — > j*q{qg — > 7*9) reactions by calculating the corresponding diagrams for massive quarks and giuons, with 
the masses being distributed according to the DQPM self energies. The obtained off-shell elementary cross sections 
then have been implemented into the PHSD transport code. 

In Fig. [2] (L.h.s.) we present PHSD results for the dilepton spectrum excess over the known hadronic sources as 
produced in In + In reactions at 158 AGeV compared to the acceptance corrected data published recently by the 
NA60 Collaboration [l6j]. We find that the spectrum at invariant masses below 1 GeV is well reproduced by the p 
meson yield, if a broadening of the meson spectral function in the medium is assumed 11 1 . Thus the result of [l2j is 
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FIG. 2: Left: Acceptance corrected mass spectra of the excess dimuons from In + In at 158 AGeV integrated over pr in 
0.2 < pr < 2.4 GeV from PHSD compared to the data of NA60 16]. Right: The inverse slope parameter T e ff of the dimuon 
yield from In+In at 158 A- GeV as a function of the dimuon invariant mass in PHSD compared to the data of the NA60 
Collaboration [l-l . fl6| . The contribution from charm decays has been subtracted from the data. 



confirmed by the PHSD calculations. On the other hand, the spectrum at M > 1 GeV is shown to be dominated by 
the partonic sources. 

It is also interesting to note that accounting for partonic dilepton sources allows to reproduce in PHSD (cf. Fig. [2] 
rhs) the intriguing finding of the NA60 Collaboration [3, [n| that the effective temperature of the dileptons (slope 
parameters) in the intermediate mass range is lower than the T e // of the dileptons from the hadronic phase. The 
softening of the transverse mass spectrum with growing invariant mass implies that the partonic channels occur 
dominantly before the collective radial flow has developed. In total we still underestimate the slope parameter T e ff 
(cf. rhs in Fig. 2) which might be due to missing partonic initial state effects or an underestimation of partonic flow 
in the initial phase of the reaction. These issues are presently being addressed and will be reported in the near future. 



IV. RHIC ENERGIES 



In 2008, the PHENIX Collaboration has presented first dilepton data irompp and Au + Au collisions at Relativistic- 
Heavy-Ion-Collider (RHIC) energies ofVi= 200 GeV [3] which show an even larger enhancement in Au+Au reactions 
(relative to pp collisions) in the invariant mass regime from 0.15 to 0.6 GeV than the data at SPS energies [? ]. We 
recall that HSD provides a reasonable description of hadron production in Au + Au collisions at y/s = 200 GeV [2(| 
such that we can directly continue with the results for e + e~ pairs which are shown in Fig. [3] for p + p collisions. We 
find that the dilepton radiation in the elementary channel is well under control in HSD, which is equivalent to PHSD 
for p + p reactions. 

We step on with the case of inclusive Au + Au collisions in comparison to the data from PHENIX [lj| (Fig [3| . 
When including the in-medium modification scenarios for the vector mesons, we achieve a sum spectrum which is 
only slightly enhanced compared to the 'free' scenario. Whereas the total yield is quite well described in the region 
of the pion Dalitz decay as well as around the w and <f> mass, HSD clearly underestimates the measured spectra in 
the regime from 0.2 to 0.6 GeV by an average factor of 3. The low mass dilepton spectra from Au + Au collisions at 
RHIC (from the PHENIX Collaboration) are clearly underestimated in the invariant mass range from 0.2 to 0.6 GeV 
in the 'collisional broadening' scenario as well as in the 'dropping mass + collisional broadening' model, i.e. when 
assuming a shift of the vector meson mass poles with the baryon density. We mention that our results for the low 
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FIG. 3: The HSD results for the mass differential dilepton spectra in case of pp collisions at y/s = 200 GeV in comparison to 
the data from PHENIX [T^. The actual PHENIX acceptance and mass resolution have been incorporated (see legend for the 
different color coding of the individual channels). Figure taken from [ill ]. 



mass dileptons are very close to the calculated spectra from van Hees and Rapp as well as Dusling and Zahed [2l[ (cf. 
the comparison in Ref. [22]). Consequently we attribute this additional low mass enhancement seen by PHENIX to 
non-hadronic sources and will proceed with the quantitative test of this assumption in a future publication. 
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FIG. 4: The HSD results for the mass differential dilepton spectra in case of inclusive Au + Au collisions at t/s = 200 
GeV in comparison to the data from PHENIX 18]. The actual PHENIX acceptance filter and mass resolution have been 
incorporated [l{|. In the upper part the results are shown for vacuum spectral functions (for p,Lu,(j>) including the channel 
decompositions (see legend for the different color coding of the individual channels). The lower part shows a comparison for 
the total e + e~ mass spectrum in case of the 'free' scenario (dashed line), the 'collisional broadening' picture (solid line) as well 
as the 'dropping mass + collisional broadening' model (dash-dotted line). Figure taken from [ll[. 
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